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Abstract 

We present first results for masses of spin-| and A baryons in lattice QCD using Fat-Link 
Irrelevant Clover (FLIC) fermions. Spin-| interpolating fields providing overlap with both spin-| 
and spin-i states are considered. In the isospin-i sector, we observe, after appropriate spin and 
parity projection, a strong signal for the J^ = | state together with a weak but discernible 
signal for the | state with a mass splitting near that observed experimentally. We also find 
good agreement between the ^ masses and earlier nucleon mass simulations with the standard 
spin-^ interpolating field. For the isospin-| A states, clear mass splittings are observed between 
the various 2 and ^ channels, with the calculated level orderings in good agreement with those 
observed empirically. 
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I. INTRODUCTION 



The level orderings in the baryon spectrum and mass splittings between excited baryon 
multiplets provide important clues to the underlying dynamics governing inter-quark forces 
and the relevant effective degrees of freedom at low energy Considerable insights into 
these and other problems of spectroscopy have been gained from QCD-inspired phenomeno- 
logical models, however, many fundamental questions about the origins of the empirical 
spectrum remain controversial |j2|]. 

The resolution of some of these issues may only be possible with the help of calculations of 
the spectrum in lattice QCD — currently the only first-principles method able to determine 
hadron properties directly from the fundamental quark and gluon theory. Recent advances 
in computational capabilities and more efficient algorithms have enabled the first dedicated 
lattice QCD simulations of the excited states of the nucleon to be undertaken 0, 0, S S 
H, 1^ • Lattice studies of excited hadrons are possible because at the current unphysically 
large quark masses and finite volumes used in the simulations, most excited states are stable. 
Contact with experiment can be made via extrapolations incorporating the nonanalytic 



behaviour of chiral effective field theory [1J|. These studies are timely as they complement 
the first results from the high precision measurements of the A^* spectrum at Jefferson Lab 

131. 



In a recent paper ^ we presented first results for the excited nucleon and spin-| hy- 
peron spectra using the Fat-Link Irrelevant Clover (FLIC) quark action |13|] with an (9(a^)- 
improved gluon action. The FLIC action minimizes the effect of renormalization of action 
improvement terms and displays excellent scaling properties Clear mass splittings are 



observed for = and | states, as well as evidence for the sensitivity to hyperfine 
splittings of odd-parity states for different interpolating fields used. On the other hand, no 
evidence is seen for overlap of three-quark interpolating fields with the Roper-like resonances 
or the lowest-lying odd-parity SU(3) singlet state, the A (1405). In this paper we extend the 
analysis of Ref. [9| to the spin-| sector, and present first results using the FLIC action, in 
both the isospin-i and | channels. 

Mass splittings between states within SU(3) quark- model multiplets provide another 
important motivation for studying higher spin baryons. Understanding the mass split- 
ting between the (1535) and (1520), for instance, or between the A| (1620) and 
A| (1700), can help identify the important mechanisms associated with the hyperfine in- 
teract ions^or shed light on the spin-orbit force, which has been a central mystery in spec- 
troscopy [ij]. In valence quark models, the degeneracy between the and can be 
broken by a tensor force associated with mixing between the iV^ and N"^ representations of 
SU(3) [^, although this generally leaves the at a higher energy than the . On 

the other hand, a spin-orbit force is necessary to split the A| and A| states. In the 
Goldstone boson exchange model [l^, both of these pairs of states are degenerate. Model- 
independent analyses in the large Nc limit have found that these mass splittings receive 
important contributions from operators that do not have a simple quark model interpreta- 
tion |16], such as those simultaneously coupling spin, isospin and orbital angular momentum. 
Of course, the coefficients of the various operators in such an analysis must be determined 
phenomenologically and guidance from lattice QCD is essential. 

Our lattice simulations are performed on the Orion computer cluster dedicated to lattice 
QCD at the CSSM, University of Adelaide. In the isospin-| sector, after applying suitable 

parity and spin projections, we present the first results for the A|^ and A| states, as well 
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as the A I . Our results for the A| are also in good agreement with earlier simulations 

A significant advance of this work is the observation of a discernible signal for the 
state, which yielded a weak signal in earlier simulations 0]. The lowest excitation of the 
ground state, namely the A^ , is found to have a mass ~ 350-400 MeV above the A| , 

with the A^ slightly heavier. The A-^ state is found to lie ~ 100-200 MeV above these, 
although the signal becomes weak at smaller quark masses. This level ordering is consistent 
with that observed in the empirical mass spectrum. 

In the spin-| nucleon sector, there is good agreement for the spin-projected and | 
states with earlier nucleon mass calculations using the standard spin-| nucleon interpo- 
lating field. Furthermore, we find a good signal for the states, with a mass difference of 

~ 300 MeV between the spin-| parity partners. The and states are approximately 
degenerate as observed experimentally. 

In Section |n] we outline the basic elements of formulating spin-| baryons on the lattice, 
including the choice of interpolating fields and projection operators. A brief prelimi nary 
report of states using the formalism developed and presented here appeared in Ref. jl8| . 
These results also supersede preliminary results reported in Refs. [19. , ,2ol | . In Section llll|, 
our results using the FLIC action on a large lattice volume at a fine lattice spacing represent 
the first quantitative analysis of these states. The conclusion and remarks about future work 
are contained in Section HVl 



II. SPIN 3/2 BARYONS ON THE LATTICE 

A. Spin 3/2 Interpolating Fields and Two-Point Functions 

In this Section the essential elements for a lattice calculation of spin-| baryon properties 
are presented. The mass of a spin-| baryon on the lattice is obtained from the two-point 
correlation function G^,^ jl7l |. 

G^,{t,p;T)=ti,^{Tg^,{t,p)} , (1) 

where 

0"J{t,p) = X^e-P- (0| T (x^(a:) x^lO)) |0) , (2) 

X 

where is a spin-| interpolating field, F is a matrix in Dirac space with a, (3 Dirac indices, 
and /X, v Lorentz indices. 

In this analysis we consider the following interpolating field operator for the isospin-|, 
spin-|, positive parity (charge -|-1) state [2l| 

^ ^a,c (^Ta^^) ^^^^u ^6^^)) ^^^^ _ 1^^^^^ ^^^c^^^ _ (3) 

All discussions of interpolating fields are carried out using the Dirac representation of the 7 
matricies. This exact isospin-| interpolating field has overlap with both spin-| and spin-| 
states and with states of both parities. The resulting correlation function will thus require 
both spin and parity projection. The quark field operators u and d act at Euclidean space- 
time point X, C is the charge conjugation matrix, a, h and c are colour labels, and the 
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superscript T denotes the transpose. The charge neutral interpolating field is obtained by 
interchanging u ^ d. This interpolating field transforms as a Rarita-Schwinger operator 
under parity transformations. That is, if the quark field operators transform as 

Vu{x)V^ = +7om(x) , 

where x = (xq, —x), and similarly for d{x), then 

and similarly for the Rarita-Schwinger operator 

Pu^ix)V^ = +7ou^{x) . (4) 

which will be discussed later. 

The computational cost of evaluating each of the Lorentz combinations in Eq. Q is 
relatively high — about 100 times that for the ground state nucleon |4]. Consequently, in 
order to maximize statistics in our analysis we consider only the leading term proportional 
to Qfj^i, in the interpolating field, 

e'-'' {u^'^ix) d\x)) ^,u\x) . (5) 

This is sufficient since we will in either case need to perform a spin-| projection. 

In order to show that the interpolating field defined in Eq. (jS)) has isospin-^, we first 
consider the standard proton interpolating field, 

X" = e'''\u^''C^,d')u\ (6) 

which we know to have isospin-|. Applying the isospin raising operator, on x^ one finds, 

= 0. 

Similarly, for the interpolating field defined in Eq. (0), one has 

= -e'^^^(M^"C757Mw')75W^ 
= 0, 

where we have used the representation independent identities C7^C^^ = — 7J, Cj^C'^ = ■jj 
and the identities which hold in the Dirac representation C"^ = = C^^ = —C with 
C = i727o and 7^ = 75- 

We note that x^ corresponding to x^ iii Eq. © is 

X^ = e^'^u^l.id'^.l.Cu^^) , (7) 
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so that 



75S'„75tr 



where the last hne is the result achieved after doing the Grassman integration over the quark 
fields with the quark fields being replaced by all possible pairwise contractions. 

In deriving the A interpolating fields, it is simplest to begin with the state containing 
only valence u quarks, namely the A"*"^. The interpolating field for the A"^"*" resonance is 
given by (2 In] . 

^(x) = e^'^ (n^«(a;) Cj, u\x)) u%x) , (9) 

which also transforms as a pseudovector under parity. The interpolating field for a A+ state 
can be similarly constructed 117^ . 



xf (x) = i= e'^^^ [2 {u^'^ix) d\x)) u\x) + [u'^ix) C^, u\x)) d\x)\ . (10) 

Interpolating fields for other decuplet baryons are obtained by appropriate substitutions of 
u, d u, d or s fields. 

To project a pure spin-| state from the correlation function G^j,^, one needs to use an 
appropriate spin-| projection operator 2^ |. 

P^'^ (P) = 9i^u - '^7^,7u - ■Pl^^Pu+P^,7ul -p) ■ (H) 

The corresponding spin-| state can be projected by applying the projection operator 

P',i' = 9.. - PT ■ (12) 

To use this operator and retain all Lorentz components, one must calculate the full 4x4 
matrix in Dirac and Lorentz space. However, to extract a mass, only one pair of Lorentz 
indices is needed, reducing the amount of calculations required by a factor of four. We 
calculate the third row of the Lorentz matrix and use the projection, 

4 

^33= 5^G3,^^^P;3, (13) 

to extract the desired spin states, s = | or |. Following spin projection, the resulting 
correlation function, Ggg, still contains positive and negative parity states. 



B. Baryon Level 

The interpolating field operators defined in Eqs. Q and (0) have overlap with both 
spin- 1 and spin-i states with positive and negative parity. The field x^i transforms as a 
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pseudovector under parity, as does the Rarita-Schwinger spinor, u^. Thus the overlap of 
with baryons can be expressed as 



'My 2 



{0\x,\N-^^ip,s)) = Xy,^J-^ u,ip,s) , (14a) 



E3/2- 
I Ms/2 



{Q\x,\N-^-ip,s)) = Xy.^J-^ 75U,ip,s) , (14b) 



-£■3/2 



{Q\Xt^\N'^^{p,s)) = (ai/2+PM + A/2+7m)W^^ 75m(p,s) , (14c) 

V -^1/2+ 

1 / Ml /2- 

{0\Xf.\N-2~{p,s)) = (ai/2-p^ + /3i/2-7M)WTr^ ^(P^^) ' (14d) 

V -^1/2- 

where the factors A5, a^, /J^ denote the couphng strengths of the interpolating field X/i 

to 

the baryon B and Eb = \/p'^ + is the energy. For the expressions in Eqs. ()14cj) and 
()14d|) . we note that the spatial components of momentum, pi, transform as a vector under 
parity and commute with 70, whereas the 7^ do not change sign under parity but anticommute 
with 7o. Hence the right-hand-side of Eq. p4c|) also transforms as a pseudovector under 
parity in accord with X/i- 

Similar expressions can also be written for x^. 



(iVt+(p,s)|xjO) = A;/2+./^^ u,{p,s) , (15a) 

\/ -&3/2+ 



'M3/2 



{N---{p,s)\x,\0) = -K/2-\l^^ Mp,sH , (15b) 



{N^^{p,s)\x^,\0) = -W^^^ m(p, 5)75(0^/2+ + A*/2+^A') ' (^^^) 

y -£^1/2+ 

{N---{p,s)\x,\0) = J^^u{P,s){al^2-P, + (3*y,^^,) . (15d) 

Note that we are assuming identical sinks and sources in these equations. In our calculations 
we use a smeared source and a point sink in which case A*, a* and (3* are no longer complex 
conjugates of A, a and /3 and are instead replaced by A, a and /?. 

We are now in a position to find the form of Eq. (0) after we insert a complete set of 
intermediate states {\B{p, s))}. The contribution to Eq. Q from each intermediate state 
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considered is given by 

{0\x>.\nI+{p,s)){nI+{p,s)\xu\0) 

^ M3/2 



+A3/2+A3/2+ Ui,{p, s)u„{p, s) 



E3/2 



- M3/2+ (7-P + M3/2+) j 1 2p^p^ p^llu-Pvl^Ji 
-A3/2+A3/2+ TTT \ 9tiu - 1:1 ^Hv - ^jj2 ^ 



E3/2+ 2M3/2+ [''^^ 3'^"" 3M3^/2+ 3M3/24 



{0\x,\N-^-iP,s)){N"(p,s)\xM 



M3/2 



A3/2-A3/2- IbU^^ip, s)uy{p, s)75 



-E3/2 



^ Y ^^3/2- (7 ■ P - ^3/2- ) j 1 2p^p^ p^7^ - Pj,7^ 

-A3/2-A3/2- TTTZ \ g^iv - -Itilu 



E3/2- 2M3/2- r'^^ 3"^"^ SMgy 3M3/2- 



{0\x,\N-^^ip,s)){m+ip, s)\xu\0) 



+ Pi/2+7m)75 — 75(ai/2+Pz. + Pi/2+ly) 



-£^1/2+ ZiWi/2+ 

Ml /2- . 7 ■ P + Ml /2- ,_ - 

= [0ll|2-P^, + P1/2-I1.) 7^ [Oiil2-Pu + ■ 

-C/-/2+ /Mi/2- 

To reduce computational expense, we consider the specific case when = z/ = 3 and in order 
to extract masses we require p = (0, 0, 0). In this case we have the simple expressions 

(0|x3|iVi+(p,.))(iVf+(p,.)|x3|0) = A3/2.A3/2+^ { ^'^'2M^/2^'^'' ) ' ^^^^^ 
{0\X3\nI-{p,s)){nI-{p,s)\x3\0) = X3/2~Xv2^l { ^'^'2M~/2^'''~ ) ' ^^^^^ 
(0|x3|iV^+(p,^))(iV^+(p,^)|X3|0) = -A/2+^i/2+7375 ^°^y" ^ 7573 

-75 70M1/2+ + M1/2+ 

= +A/2./?V2+ , (16c) 

{0\X3\N'^-{p,s)){N'^~ip,s)\xM = ^1/2-^,2-13 ^'^"':,'^ 13 

7^ 70M1/2- - Ml/2- 
= +/5i/2-/3v2- . (16d) 

Therefore, in an analogous procedure to that used in Ref. P, where a fixed boundary 
condition is used in the time direction, positive and negative parity states are obtained by 
taking the trace of the spin-projected correlation function, G33, in Eq. (fTSj) with the operator 

r = r±, 

G^f = tr,p{r±G^3} , (17) 

where 

r± = ^(i±74). (18) 
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The positive parity states propagate in the (1,1) and (2,2) elements of the Dirac matrix, 
while negative parity states propagate in the (3,3) and (4,4) elements for both spin-| and 
spin-| projected states. A similar treatment has been carried out for the A interpolating 
fields but is not reproduced here for brevity. 

We use an improved unbiased estimator obtained by summing both U and U* configu- 
rations which occur with equal weight. From the discussion given in Section Va of Ref j9|, 
is purely real if fi and u are both spatial indicies or both temporal indicies, otherwise 
Cfj^ is purely imaginary. 



III. RESULTS 

The analysis is based on a sample of 392 configurations. For the gauge fields, a mean- 
field improved plaquette plus rectangle action is used. The simulations are performed on a 
16^ X 32 lattice at (3 = 4.60, which corresponds to a lattice spacing of a = 0.122(2) fm set by a 
string tension analysis incorporating the lattice coulomb potential ji^ with y/a = 440 MeV. 
For the quark fields, a Fat-Link Irrelevant Clover (FLIC) action is implemented. The 
use of fat links 2^1 in the irrelevant operators of the fermion action removes the need to 



fine tune the clover coefficient to remove 0{a) artifacts. The use of fat links also allows us 
to employ a highly improved definition of Ff^,^ jisl . |25| leaving errors of 0{a^) and where 
errors of 0{g^) are suppressed via fat links. Mean- field improvement of the tree-level clover 
coefficient with fat links represents a small correction and proves to be adequate [l3| . 

The fattening, or smearing of the lattice links with their nearest neighbours, reduces the 
problem of exceptional configurations, and minimizes the effect of renormalization on the 
action improvement terms. By smearing only the irrelevant, higher dimensional terms in 
the action, and leaving the relevant dimension-four operators untouched, we retain short 
distance quark and gluon interactions at the scale of the cutoff. Simulations are performed 
with n = 4 smearing sweeps and a smearing fraction a = 0.7 jl3l] . 

A fixed boundary condition in the time direction is used for the fermions by setting 
Ut{x,Nt) = V a; in the hopping terms of the fermion action, with periodic boundary 
conditions imposed in the spatial directions. Gauge- invariant gaussian smearing j2^ in 
the spatial dimensions is applied at the source to increase the overlap of the interpolating 
operators with the ground states. The source-smearing technique j26| starts with a point 
source, 

iPoa{S,t) = 6^'6o.^6g^g,6t,to (19) 
for source colour c, Dirac 7, position xq = (1, 1, 1) and time to = 3 and proceeds via the 



iterative scheme, 



where 



i)i{x, t) = ^ F{x, f ) ipi-i{x', t) 



Fix, X) = -^Y^^ (^Sx,x' + -^Y1 [^Z^'^^' ^) + ^li^ - ^) ^x',x-jl\ 

Repeating the procedure N times gives the resulting fermion source 

^|J^{x, t)=J2 F^'i^^ M^'^ t) ■ (20) 
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TABLE I: 


Masses of the vr 


, iVi^ and A^l^, 


for several values 


of K obtained 


from the spin-| 


interpolatin 


g field, for the FLIC action with 4 


sweeps of smearin; 


5 at a = 0.7. 


Here the value of 


^cr IS — 


0.1300. A string 


tension analysis provides a = 0.122(2 


) fm for -y/o" = 


440 MeV. 


K 


m^r a 




Mj^i a 

9. 




M«r3 a 

'?. 


0.1260 


0.5767(11) 


1.102(8) 


1.412(13) 


1.628(34) 


1.410(16) 


0.1266 


0.5305(12) 


1.043(9) 


1.369(14) 


1.577(38) 


1.365(19) 


0.1273 


0.4712(15) 


0.970(13) 


1.317(17) 


1.510(44) 


1.312(24) 


0.1279 


0.4164(15) 


0.905(18) 


1.271(21) 


1.440(53) 


1.264(32) 


0.1286 


0.3421(18) 


0.829(32) 


1.220(31) 


1.329(74) 


1.206(49) 



The parameters and a govern the size and shape of the smearing function. We simulate 
with = 20 and a = 6. The propagator, 5, is obtained from the smeared source by solving 

St, = C , (21) 

for each colour, Dirac source c, 7 respectively of Eq. (fT^ via the BiStabilised Conjugate 
Gradient algorithm ^7|. 

In the analysis we use five values of as indicated in Table UJ Extrapolation to = 
gives = 0.1300. Figure [T] shows the effective mass plot for the state for the five 
K values used as a function of Euclidean time obtained after performing spin and parity 
projections on the correlation functions calculated using the interpolating field in Eq. 
We find a good signal for this state up until time slice 13 after which the signal is lost in 
noise. The effective mass for this state exhibits good plateau behaviour and a good value of 
the covariance-matrix based x^/^df is obtained when one fits in the time fitting window of 
t = 10-13 (recall, the source is at t = 3). Typically, one finds x^/^df ~ 1 and x^/^df < 1-5 
throughout. After performing spin and parity projections to extract the A^f"*" state from the 
interpolating field in Eq. (0), one finds the effective mass plot to be a little noisier, as shown 
in Fig. 121 There is, however, sufficient information here to extract a mass, and a good value 
of x^/^BF is obtained when one fits in the small time fitting window of t = 9-11. 

The interpolating field defined in Eq. (0) also has overlap with spin-i states of both 
parities. After performing a spin-^ projection on the correlation functions, we isolate the 
A^i^ and states via parity projection and plot the effective masses in Figs. El and El 
respectively. The A^^"*" state suffers contamination from excited states as seen by the long 
Euclidean time evolution required to reach plateau in Fig. El A good value of x^/^df is 
obtained as long as we fit after time slice 12. For this reason, we use time slices 13-16 to 
obtain a mass for the A^l"*^ state. However, for the state, a plateau is seen at early 
Euclidean times and a good value of x^/^df is obtained on time slices 8-11. 

The extracted masses of the A^|^ and states are given in Table U and are displayed 
in Fig. El as a function of m^. Earlier results for the A^| states using the standard spin- 
2 interpolating field 0, El] from Eq. © are also shown with open symbols in Fig. El for 
reference. It is encouraging to note the agreement between the spin-projected | states 
obtained from the spin-| interpolating field in Eq. (0) and the earlier ^ results from the 
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FIG. 1: Effective mass plot for the state using the FLIC action, from 392 configurations. The 
five sets of points correspond to the k values listed in Table with k increasing from top down. 
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FIG. 2: As in Fig.d but for the iV|"^ state. 
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FIG. 3: As in Fig.Ul but for the iVi^ state. 



1.8 
1.6 

1.4 h 

b 1.2 
1.0 
0.8 



♦ - 





i' L 



3 4 5 6 7 



8 9 10 11 12 13 14 
t 



FIG. 4: As in Fig. [U but for the N\ state. 
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FIG. 5: Masses of the spin projected (filled triangles), -/V|^ (filled inverted triangles), iV^"*" 

(filled circles), and (filled squares) states. For comparison, previous results from the direct 

calculation of the N}^^ (open circles) and N}^ (open squares) from Ref. 0| are also shown. The 
empirical values of the masses of the N\'^ (939), (1535), N^' (1520) and iVf (1720) are 

shown on the left-hand-side at the physical pion mass. 



same gauge field configurations. To study this agreement more accurately, we consider the 
ratio of effective masses obtained for each jackknife subensemble. This provides us with a 
correlated ratio and we find the ratio to be one as at the one standard deviation level. We 
also observe that the state has approximately the same mass as the spin-projected 

iVi state which is consistent with the experimentally observed masses. To study this 
mass difference more accurately, we again calculate the correlated ratio of effective masses 
obtained after appropriate spin and parity projections. This ratio is found to be one within 
one standard deviation. The results for the state in Fig. El indicate a clear mass splitting 

between the A^|^ and states obtained from the spin-| interpolating field, with a mass 
difference around 300 MeV. This is slightly larger than the experimentally observed mass 
difference of 200 MeV. 

Turning now to the isospin-| sector, the effective mass plot for the Al"*" state using the 
interpolating field given in Eq. is shown in Figure El for the five k values used. An 
excellent signal is clearly visible, and a good value of the covariance-matrix based x^/^bf 
is obtained by fitting time slices t = 11-14 following the source at t = 3. For the effective 
mass of the negative parity A| , shown in Fig. [71 the signal is quite good up to time slice 
11-12, but is lost in noise after time slice 12. Time slices t = 9-12 provide a fitting window 
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FIG. 6: Effective mass plot for the state using the FLIC action with 4 sweeps of smearing 

at a = 0.7 from 392 configurations. The five sets of points correspond to the k values listed in 
Table ^ with k increasing from top down. 



TABLE II: As in Tabled but for the corresponding 



A 2 and A 2^ 



masses. 





M 3^ a 


M 3 a 


M i.a 


M 1 a 


0.1260 


1.198(8) 


1.469(15) 


1.643(109) 


1.476(34) 


0.1266 


1.153(9) 


1.429(17) 


1.604(107) 


1.432(41) 


0.1273 


1.101(12) 


1.385(21) 


1.561(106) 


1.387(54) 


0.1279 


1.057(15) 


1.353(27) 


1.530(109) 


1.351(76) 


0.1286 


1.006(22) 


1.331(43) 


1.502(119) 


1.301(126) 



with an acceptable value of x^/^df- 

The results for the Al"*^ and A| masses are shown in Fig. |Hlas a function of m^. The 
trend of the A^~^ data points with decreasing is clearly towards the A(1232), although 



some nonlinearity with is expected near the chiral limit |28l. I29| . The mass of the A| 



lies some 500 MeV above that of its parity partner, although with somewhat larger errors, 
as expected from the effective mass plots in Figs. El and [3 

After performing a spin projection to extract the A^ states, a discernible, but noisy, 
signal is detected. This indicates that the interpolating field in Eq. Q has only a small 
overlap with spin-^ states. However, with 392 configurations we are able to extract a mass 
for the spin- 1 states at early times, shown in Fig. |H1 Here we see the larger error bars 



13 



\ r 



1.8 
1.6 
^ 1.4 
1.2 
1.0 
0.8 



J L 



/ i 



3 4 5 6 



7 8 9 10 11 12 13 14 
t 



FIG. 7: As in Fig. H but for the A| state. 



associated with the states. The lowest excitation of the ground state, namely the , 

has a mass ~ 350-400 MeV above the A| , with the A| possibly appearing heavier. The 

Ai"*^ state is found to lie ~ 100-200 MeV above these, although the signal becomes weak at 
smaller quark masses. This level ordering is consistent with that observed in the empirical 
mass spectrum, which is also shown in the figure. 

The iV| and A| states will decay to iVvr in S'-wave even in the quenched approximation 
[3^ . For all quark masses considered here, with the possible exception of the lightest quark, 
this decay channel is closed for the nucleon. While there may be some spectral strength in 
the decay mode, we are unable to separate it from the resonant spectral strength. 

The iV|^ and A^"*" states will decay to iVvr in P wave, while iV| and A| states will 
decay to iVvr in D-wave. Since the decay products of each of these states must then have 
equal and opposite momentum and energy given by 

these states are stable in our calculations. 



IV. CONCLUSION 

We have presented the first results for the spectrum of spin-| baryons in the isospin-| 
and I channels, using a novel Fat Link Irrelevant Clover (FLIC) quark action and an O(a^) 

si 1 it 

improved gauge action. Clear signals are obtained for both the spin-projected and N ^ 
states from a spin-| interpolating field. In particular, the states are in good agreement 
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FIG. 8: Masses of the spin-projected and A^^ resonances. The empirical values of the masses 
of the A3/2+(l232), A|~(1700), Ai~(1620) and Ai+(1910) are shown on the left-hand-side at the 
physical pion mass. 



with earlier simulations of the nucleon mass and its parity partner using the standard spin- 
I interpolating field. We find the state to he at a similar energy level to the , 
consistent with experiment. We also find a mass difference of ~ 300 MeV between the 
spin- 1 , isospin-| parity partners, slightly larger than the experimentally observed difference 
of 200 MeV. 

For isospin-| baryons, good agreement is found with earlier calculations for the A ground 
state, and clear mass splittings between the ground state and its parity partner arc observed 
after suitable spin and parity projections. Wc obtain a signal for the A| states and the 
level ordering is consistent with that observed in the empirical mass spectrum. 

It will also be important in future work to consider the excited states in each channel, 
in particular the lowest "Roper-like" excitation of the A(1232) ground state. Although this 
will be more challenging, it may reveal further insights about the origin of the inter-quark 
forces and the nature of the confining potential. 
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